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Thermalization in strong interactions 



1. Introduction 

Quantum chromodynamics, through colour confinement, restricts quarks and gluons to move 
in a limited region of space. In a black hole, matter and light are confined by gravitation to remain 
within a restricted region of space. This similarity was noted quite soon after the advent of QCD, 
and it was suggested that hadrons were the analogue of black holes in strong interaction physics 
[jj], ||, 0|. In the case of gravitational black holes, the isolation of the system is not quite absolute: in 
the strong field at the outer edge, quantum excitation can lead to the emission of Hawking radiation 
into the physical vacuum [Qj. Since no information transfer between the inside and the outside of 
the black hole is allowed, this radiation must give equal a priori weights to all possible states on the 
outside, and hence it is thermal at the point of formation. In this report, based on joint work with 
Paolo Castorina and Dmitri Kharzeev [|5j|, we want to show that high energy collisions in strong 
interaction physics produce a self-similar cascade of "white holes" (colourless from the outside, 
but coloured inside), and that the Hawking radiation arising at the confinement horizon of these 
can provide the thermal behaviour observed quite universally in all soft hadron production |7|] . 

We begin by recalling briefly the relevant basic aspects of black hole physics and indicate 
some first relations to strong interactions. We then turn to high energy e + e~ annihilation as the 
simplest hadron production process and study the effect of quark-antiquark pair excitation and 
string breaking. Next we review some features of charged and of rotating black holes; these allow 
us to generalize our scenario to hadronic collisions at finite baryon density as well as to non-central 
interactions. Finally, we comment on stochastic vs. kinetic thermalization. 



2. Black Holes and Hawking Radiation 

A black hole is the final stage of a neutron star after gravitational collapse It has a mass 
M concentrated in such a small volume that the resulting gravitational field confines all matter 
and even photons to remain inside the event horizon R of the system: no causal connection with 
the outside is possible. As a consequence, black holes have three (and only three) observable 
properties: mass M, charge Q and angular momentum J. This section will address mainly black 
holes with Q = J = 0; we shall come back to the more general properties in section 4. 

Classically, a black hole would persist forever and remain forever invisible. On a quantum 
level, however, its constituents (photons, leptons and hadrons) have a non-vanishing chance to 
escape by tunnelling through the barrier presented by the event horizon. The resulting Hawking 
radiation [Q] cannot convey any information about the internal state of the black hole; it must 
be therefore be thermal, and it was shown that for a non-rotating black hole of vanishing charge 
(denoted as Schwarzschild black hole), the corresponding radiation temperature is 

Tbh = ; (2-1) 

&71GM 

where G is the gravitational constant. This temperature is inversely proportional to the mass of the 
black hole, and since the radiation reduces the mass, the radiation temperature will increase with 
time. For black holes of stellar size, however, one finds Tbh <2x 10~ 8 °K. This is many orders of 
magnitude below the 2.7 °K cosmic microwave background, and hence not detectable. 



2 



Thermalization in strong interactions 



In general relativity, the event horizon arises as consequence of the Schwarzschild metric and 
its generalizations to Q ^ 0,7 7^ 0. The occurrence and role of the event horizon for thermal 
radiation was subsequently generalized by Unruh [^]. A system undergoing uniform acceleration 
a relative to a stationary observer eventually reaches a classical turning point and thus encounters 



an event horizon. Let us recall the resulting hyperbolic motion (JlOJ. A point mass m subject to a 
constant force F satisfies the equation of motion 



dty/Y 



(2.2) 



where v(?) = dx/dt is the velocity, normalized to the speed of light c = 1. This equation is solved 
by the parametric form 

jc = -coshaT ? = -sinhaT, (2.3) 



where a = F/m denotes the acceleration in the instantaneous rest frame of m, and T the proper 
time, with dx = y/T— v 2 dt. 
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Figure 1: Hyperbolic motion 

The resulting world line is shown in Fig. [jj It corresponds to the mass m coming from x = 00 
at t = —00 at with a velocity arbitrarily close to that of light, decelerating uniformly until it comes 
to rest at the classsical turning point x# = — (1/a). It then accelerates again and returns to x = 00 
at t = 00, approaching the speed of light. For given a, xh is thus the event horizon beyond which 
m classically cannot pass. The only signal an observer can detect of the passage of m is thermal 
quantum radiation of temperature 



11, 12] 



a 

2n' 



In the case of gravity, we have the force 



ma 



G- 



mM 



(2.4) 



(2.5) 



on a probe of mass m. From the dependence of the gravitational metric on black hole mass and 
radius, we obtain 



R e = 2GM 



(2.6) 



for the Schwarzschild radius of the black hole, leading to a = 1/(4 GM) for the acceleration at the 



event horizon and thus back to eq. (2.1). We also see why the temperature of Hawking radiation 
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decreases with increasing black hole mass: since R g grows linearly with M, increasing the black 
hole mass decreases its energy density and hence the Hawking temperature Tbh- 

It is instructive to consider the Schwarzschild radius of a typical hadron, assuming a mass 
m ~ 1 GeV: 

R h g ad ~ 1.3 x 1(T 38 GeV" 1 ~ 2.7 x 1(T 39 fm. (2.7) 

To become a gravitational black hole, the mass of the hadron would thus have to be compressed 
into a volume more than 10 100 times smaller than its actual volume, with a radius of about 1 fm. On 
the other hand, if instead we increase the interaction strength from gravitation to strong interaction 
[ffl], we gain in the resulting "strong" Schwarzschild radius R^ ad a factor 

(X " (2.8) 



Gm 2 ' 

where a s is the dimensionless strong coupling and Gm 2 the corresponding dimensionless gravita- 
tional coupling for the case in question. This leads to 

R h s ad ~ ^ (2.9) 
m 

which for a s ~ 2.5 gives R h ° d = 1 fm. In other words, the confinement radius of a hadron is about 
the size of its "strong" Schwarzschild radius, so that we could consider quark confinement as the 
strong interaction version of the gravitational confinement in black holes [Q,^]. 

The event horizon of a black hole arises from the deformation of the Minkowski metric due 
to the strong gravitational force. Such a geometrization of force effects has also been employed in 



non-linear electrodynamics, where it can lead to similar spatial constraints []13[]. A corresponding 
extension to QCD will be discussed in ; we shall here concentrate on the more phenomenological 
aspects. 

3. Pair Production and String Breaking 

In the previous section, we had considered constituents of a black hole undergoing accelerated 
motion in classical space-time. In this section, we shall first address the modifications which arise 
if the underlying space-time manifold is specified by quantum field theory, so that in the presence 
of a strong field the vacuum becomes unstable under pair production. Next we turn to the specific 
additional features which come in when the basic constituents are subject to colour confinement 
and can only exist in colour neutral bound states. 

As basic starting point, we consider two-jet e + e~ annihilation at cms energy y/s, 

e + e~ — > 7 — > qq — > hadrons. (3.1) 

The initially produced qq pair flies apart, subject to the constant confining force given by the string 
tension a; this results in hyperbolic motion [14, 15, [If]] of the type discussed in the previous 



section. At t = 0, the q and q separate with an initial velocity vq = p/ \J p 2 + m 2 , where p ~ yfs/2 
is the momentum of the primary constituents in the overall cms and m the effective quark mass. 



We now have to solve Eq. ( |2.2| ) with this situation as boundary condition; the force 

F = a, (3.2) 
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is given by the string tension a binding the qq system. The solution is 



x=[\- yjl-vot + t 2 } 
with x = x/xq and f = t/xo; here the scale factor 



m 1 



x 



1 

- 7 



1 -v 2 a 
1 v 



(3.3) 



(3-4) 



is the inverse of the acceleration a measured in the overall cms. The velocity becomes 

dx (vo/2)— t 



v(t) 



it vanishes for 



thus defining 



x(f 



2 2 a 



^y(i-)/ 1 -^/ 4 ))^ 



2a 



(3-5) 



(3-6) 



(3.7) 



as classical turning point and hence as the classical event horizon measured in the overall cms (see 
Fig- 1). 




Figure 2: Classical and quantum horizons in qq separation 

Eq. (|3/7|) allows the q and the q to separate arbitrarily far, provided the pair has enough initial 
energy; this clearly violates colour confinement. Our mistake was to consider the qq system in a 
classical vacuum; in quantum field theory, the vacuum itself contains virtual qq pairs, and hence 
it is not possible to increase the potential energy of a given qq state beyond the threshold value 
necessary to bring such a qq pair on-shell. In QED, the corresponding phenomenon was addressed 



by Schwinger [17], who showed that in the presence of a constant electric field of strength $ the 
probability of producing an electron-positron pair is given by 



P{M,S) ~ exp{-7rm 2 /W}, 



(3.8) 



with m e denoting the electron mass and e denoting the electric charge. In QCD, we expect a similar 
effect when the string tension exceeds the pair production limit, i.e., when 



a x > 2m 



(3-9) 
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where m specifies the effective quark mass. Beyond this point, any further stretching of the string 
is expected to produce a qq pair with the probability 

P(M,o)~exp{-nm 2 /o}, (3.10) 

with the string tension a replacing the electric field strength eS. This string breaking acts like 
a quantum event horizon x q = 2 m/a, which becomes operative long before the classical turning 
point is ever reached (see Fig. ||). Moreover, the resulting allowed separation distance for our qq 
pair, the colour confinement radius xq, does not depend on the initial energy of the primary quarks. 

There are some important differences between QCD and QED. In case of the latter, the initial 
electric charges which lead to the field £ can exist independently in the physical vacuum, and the 
produced pair can be simply ionized into an e + and an e~ . In contrast, neither the primary quark 
nor the constituents of the qq pair have an independent existence, so that in string breaking colour 
neutrality must be preserved. As a result, the Hawking radiation in QCD must consist of qq pairs, 
and these can be produced in an infinite number of different excitation states of increasing mass 
and degeneracy. Moreover, the qq pair spectrum is itself determined by the strength a of the field, 



in contrast to the exponent m e jS in eq. (|3.8|), where the value of $ has no relation to the electron 
mass m e . 

Hadron production in e + e~ annihilation is believed to proceed in the form of a self-similar 



cascade [|18|, |19p. Initially, we have the separating primary qq pair, 

y^[qq] (3.11) 

where the square brackets indicate colour neutrality. When the energy of the resulting color flux 
tube becomes large enough, a further pair qiq\ is excited from the vacuum, 



(3.12) 



Although the new pair is at rest in the overall cms, each of its constituents has a transverse mo- 
mentum kj determined, through the uncertainty relation, by the transverse dimension rj of the flux 
tube. The slow q\ now screens the fast primary q from its original partner q, with an analoguous 
effect for the q\ and the primary antiquark. To estimate the qq separation distance at the point of 
pair production, we recall that the basic thickness of the flux tube connecting the qq pair is in string 
theory given by [ |20| ] 

r T = J—; (3.13) 

V 7lO 

higher excitations lead to a greater thickness and eventually to a divergence (the "roughening" 
transition). From the uncertainty relation we then have 

(3.14) 




With this transverse energy is included in eq. (p.9|), we obtain for the pair production separation xq 
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with a = 0.2 GeV 2 and m 2 <C a. 

Once the new pair is present, we have a colour-neutral system qq\q\q; but since there is a 
sequence of connecting string potentials qq\, q\q\ and q\q, the primary string is not yet broken. To 
achieve that, the binding of the new pair has to be overcome, i.e., the q\ has to tunnel through the 
barrier of the confining potential provided by q\, and vice versa. Now the q excerts a longitudinal 
force on the q\, the q on the q\, resulting in a longitudinal acceleration and ordering of q\ and q\. 
When (see Fig. |3j> 



ax{q\q\) = 2y m 2 + kj,, (3.16) 
the q\ reaches its q\q\ horizon; on the other hand, when 



Ox(qqi) =2y / m 2 +4, (3.17) 

the new flux tube qq\ reaches the energy needed to produce a further pair <72<?2- The q~2 screens the 
primary q from the q\ and forms a new flux tube qq~2- At this point, the original string is broken, and 
the remaining pair q~\q2 form a colour neutral bound state which is emitted as Hawking radiation in 
the form of hadrons, with the relative weights of the different states governed by the corresponding 
Unruh temperature. The resulting pattern is schematically illustrated in Fig. 0. 




y 

Figure 3: String breaking through qq pair production 

To determine the temperature of the hadronic Hawking radiation, we return to the original 
pair excitation process. To produce a quark of momentum kj, we have to bring it on-shell and 
change its velocity from zero tov = kj/ (m 2 + k\ ) x l 2 ~ 1 . This has to be achieved in the time of the 
fluctuation determined by the virtuality of the pair, At = 1 /AE ~ 1 jlkj. The resulting acceleration 
thus becomes 

Av / 

a= — ~2k T ~ V27TCJ ~ 1.1 GeV, (3.18) 
At 

which leads to 

(3.19) 

for the hadronic Unruh temperature. It governs the momentum distribution and the relative species 
abundances of the emitted hadrons. 

A given step in the evolution of the hadronization cascade of a primary quark or antiquark 
produced in e + e~ annihilation thus involves several distinct phenomena. The color field created 
by the separating q and q produces a further pair q\q\ and then provides an acceleration of the 
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qi, increasing its longitudinal momentum. When it reaches the q\q~\ confinement horizon, still 
another pair q2q~2 is excited; the state q~\q2 is emitted as a hadron, the q~2 forms together with the 
primary q a new flux tube. This pattern thus step by step increases the longitudinal momentum 
of the "accompanying" qi as well as of the emitted hadron. This, together with the energy of the 
produced pairs, causes a corresponding deceleration of the primary quarks q and q, in order to 
maintain overall energy conservation. In Fig. ||, we show the world lines given by the acceleration 
q~i — ► q~i + \ {qi — > qi+\) and that of formation threshold of the hadrons qiqi+x and the corresponding 
antiparticles. 



\\\ hadron i 


^ator^^^^^^^ 


quark ac( 


eleration 

X 



Figure 4: Quark acceleration and hadronization world lines 

The energy loss and deceleration of the primary quark q in this self-similar cascade, together 
with the acceleration of the accompanying partner qi from the successive pairs brings q and qi 
closer and closer to each other in momentum, from an initial separation qq\ of y/s/2, until they 
finally are combined into a hadron and the cascade is ended. The resulting pattern is shown in Fig. 

B 




Figure 5: Hadronization in e + e annihilation 



The number of emitted hadrons, the multiplicity v(s), follows quite naturally from the picture 
presented here. The classical string length, in the absence of quantum pair formation, is given 



by the classical turning point determined in eq. (3.7). The thickness of a flux tube of such an 
"overstretched" string is known [EO] ; instead of eq. (|3.13[) we have 



2 * 1 
no ^ Q 2k+l 



no 



In 2K, 



(3.20) 



where K is the string length in units of an intrinsic string vibration measure. From eq. ( |3.7| ) we thus 
get 



no 



In- 



(3.21) 



for the flux tube thickness in the case of the classical string length. Because of pair production, this 
string breaks whenever it is stretched to the length x q given in eq. (3.15); its thickness rj at this 
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point is given by eq. (3.13). The multiplicity can thus be estimated by the ratio of the corresponding 
classical to quantum transverse flux tube areas, 



vis) ~ -i- ~ 

4 



ln 



(3.22) 



and is found to grow logarithmically with the e + e~ annihilation energy, as is in fact observed 
experimentally. 

Up to now, we have considered hadron production in e + e~ annihilation, in which the virtual 
photon produces a confined coloured qq pair as a "white hole". Turning now to hadron-hadron 
collisions, we note that here two incident white holes combine to form a new system of the same 
kind, as schematically illustrated in Fig. |6| Again the resulting string or strong colour field produces 
a sequence of qq pairs of increasing cms momentum, leading to the well-known multiperipheral 
hadroproduction cascade shown in Fig. [7[ For the further extension to nuclear collisions in the 



regime of parton saturation, see [ |21| , |22| |23j]. 




Y h h 

(a) (b) 
Figure 6: "White hole" structure in e + e~ annihilation (a) and hadronic collisions (b) 



hadrons 




h h 
Figure 7: Hadronization in hadron-hadron collisions 



In the case of hadronic, and even more so for heavy ion interactions, two new elements enter. 
The resulting system can now have an overall baryon number, from B=2 in pp to B = 400 or more 
in heavy ion collisions. To take that into account, we need to consider the counterpart of charged 
black holes. Furthermore, in heavy ion collisions the resulting hadron production can be studied 
as function of centrality, and peripheral collisions will lead to an interaction region with an overall 
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angular momentum. Hence we also need to consider rotating black holes. In the next section, we 
therefore summarize the relevant feature of black holes with Q^O, 7/0. 



4. Charged and Rotating Black Holes 

As mentioned, for an outside observer the only characteristics of a black hole are its mass 
M, its charge Q, and its angular momentum /. Hence any further observables, such as the event 
horizon or the Hawking temperature, must be expressable in terms of these three quantities. 

We had noted that the event horizon of a black hole is due to the strong gravitational attrac- 
tion, which leads to a diverging Schwarzschild metric at a certain value of the spatial extension R. 
Specifically, the invariant space-time length element ds 2 is at the equator given by 

ds 2 = (1 - 2GM/R) dt 2 - — 2 1 gm/r dr\ (4.1) 

with r and t for flat space and time coordinates; it is seen to diverge at the Schwarzschild radius 
Rs = 2GM. If the black hole has a net electric charge Q, the resulting Coulomb repulsion will 
oppose and hence weaken the gravitational attraction. As a result, the corresponding form (denoted 
as Reissner-Nordstrom metric) becomes 

ds 2 = (1 - 2GM/R + GQ/R 2 ) dt 2 — dr 2 . (4.2) 

v 1 **' ' 1-2GM/R + GQ/R 2 

For this, the divergence leads to the smaller Reissner-Nordstrom radius 



R RN = GM{\ + ^\-Q 2 /GM 2 ), (4.3) 
which reduces to the Schwarzschild radius Rs for £2 = 0. The temperature of the Hawking radiation 



now becomes Qq, 24 ] 



T BH (M,Q) = T M M,0) { 2 V^Zg7g! 1 ; (4 .4) 
V ; V ' \(\ + ^/l-Q2/GM 2 ) 2 f 

its functional form is illustrated in Fig. || We note that with increasing charge, the Coulomb repul- 
sion weakens the gravitational field at the event horizon and hence decreases the temperature of the 
corresponding quantum excitations. As Q 2 — > GM 2 , the gravitational force is fully compensated 
and the black hole is gone. 

In a similar way, the effect of the angular momentum of a rotating black hole can be incor- 
porated. It is now the centripetal force which counteracts the gravitational attraction and hence 
reduces its strength. The resulting Kerr metric must take into account that in this case the rotational 
symmetry is reduced to an axial symmetry, and with 6 denoting the angle relative to the polar axis 
6 = 0, it is (at fixed longitude) given by 

, / 2GMR \ , R 2 + a 2 cos 2 d , , , , „. 



ds 2 =l i " )dt 2 - \"V„ \ dr 2 -(R 2 + a 2 cos 2 d)dd 2 . (4.5) 

1 R 2 + a 2 cos 2 dJ R 2 — 2GMR + a 2 v ; 
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Figure 8: Radiation temperature for a charged black hole 



The angular momentum of the black hole is here specified by the parameter a = J /At; for a = 0, 
we again recover the Schwarzschild case. The general situation is now somewhat more complex, 
since eq. ( |Q| ) leads to two different divergence points. The solution 



R K = GM{\ + J\-a 2 /{GM) 2 



(4.6) 



defines the actual event horizon, corresponding to absolute confinement. But the resulting black 
hole is now embedded in a larger ellipsoid 



R E = GAt ( 1 + J 1 - [a 2 / (GM) 2 ] cos 2 6 ) , 



(4.7) 



as illustrated in Fig. The two surfaces touch at the poles, and the region between them is denoted 
as the ergosphere. Unlike the black hole proper, communication between the ergosphere and the 
outside world is possible. Any object in the ergosphere will, however, suffer from the rotational 
drag of the rotating black hole and thereby gain momentum. We shall return to this shortly; first, 
however, we note that the temperature of the Hawking radiation from a rotating black hole becomes 



T BH (M,J) = T BH (M,0) 



2y/\-a 2 /(GM) 2 



(1 + y 7 ! —a 2 /(GM) 2 



(4.8) 



For a non-rotating black hole, with a = 0, this also reduces to the Hawking temperature for the 
Schwarzschild case. 



ergosphere 




( V" 


,>'"") R A 




black hole 



Figure 9: Geometry of a rotating black hole 

To illustrate the effect of the ergosphere, imagine radiation from a Schwarzschild black hole 
emitted radially outward from the event horizon. In the case of a Kerr black hole, such an emission 
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is possible only along the polar axis; for all other values of 9, the momentum of the emitted radia- 
tion (even light) will increase due to the rotational drag in the ergosphere. This effect ceases only 
once the radiation leaves the ergosphere. Since the amount of drag depends on 8, the momentum 
of the radiation emitted from a rotating black hole, as measured at large distances, will depend on 
the latitude at which it is emitted and increase from pole to equator. 

Finally, for completeness, we note that for black holes with both spin and charge (denoted as 
Kerr-Newman), the event horizon is given by 



R KN = GM (1 + y/ 1 - [Q 2 / GM 2 ] - [a 2 /(GM) 2 ]), (4.9) 
and the radiation temperature becomes [||, |4]] 

T BH {M,Q,J) = W«,0,0) ( W-( Gg ^)/( G Ay 1 

\(l + V / l-(Ge ! +« ! )/(GM) ! ) 2 +« 2 /(GM) 2 J 

The decrease of Tbh for Q ^ 0, J ^ expresses the fact that both the Coulomb repulsion and the 
rotational force counteract the gravitational attraction, and if they win, the black hole is dissolved. 



5. Vacuum Pressure and Baryon Density 

We now want to consider the extension of charged black hole physics to colour confinement 
in the case of systems with a net baryon number. In eq. ( [4~4| ) we had seen that the reduction 
of the gravitational attraction by Coulomb repulsion in a charged black hole in turn reduces the 
temperature of Hawking radiation. The crucial quantity here is the ratio Q 2 /GM 2 of the overall 
Coulomb energy, Q 2 /R, to the overall gravitational energy, GM 2 /R. Equivalently, Q 2 /GM 2 = 
Pq/Pg measures the ratio of inward gravitational pressure Pq at the event horizon to the repulsive 
outward Coulomb pressure Pq. 

In QCD, we have a "white" hole containing coloured quarks, confined by chromodynamic 
forces or, equivalently, by the pressure B of the physical vacuum. If the system has a non- vanishing 
overall baryon number, there will be a Fermi repulsion between the corresponding quarks, and this 
repulsion will provide a pressure P(ju) acting against B, with }X denoting the corresponding quark 
baryochemical potential. We thus expect a similar reduction of the hadronization temperature as 
function of ;U. To quantify this, we have to obtain the reduction of the chromodynamic force field, 
such as the string tension a, due to baryonic repulsion. As a first estimate, we start from the 
simplest picture of colour confinement and consider is an ideal gas of massless quarks and gluons, 
held together by the vacuum pressure B. At fixed temperature T and quark baryochemical potential 
/I, the overall pressure of such a system is given by 



V 90 



7 

d b + -df 



where db and df denote the degrees of freedom of gluons and quarks, respectively; for colour SU (3) 
and two quark flavours, di, = 16 and df = 12. When the inward vacuum pressure just balances the 
combined outward kinetic motion of quarks and gluons and the fermi repulsion of the quarks, we 
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have P = 0. For jx = 0, this gives 



90 

37n 2 



1/4 



(5.2) 



On the other hand, for T = 0, when the Fermi repulsion of the quarks alone balances the vacuum 
pressure, ji. becomes 



Mo = (2tt 2 B) 1/4 = 7i (^Y'to * 3 T . 



(5.3) 



With a hadronization temperature 7b ~ 0.175 GeV, this leads to a baryochemical potential jlj, = 
3jJ. ~ 1.6 GeV. In the intermediate region, where both T and /I are finite, we want to compare the 
effect of the Fermi repulsion to the vacuum pressure through the Hawking-Unruh form, i.e., we 
replace Q 2 /M 2 in eq. ( p~4| ) by jll/jUo, giving 



r(M)/zb 



(l + Vl-(M/Mo) 4 ) 2 



(5.4) 



The resulting behaviour of T(fi) is shown in Fig. 10 



1.0 



0.5 



T(u)/T(0) 




1 


^0 



0.5 1.0 

Figure 10: Hadronization temperature as function of the baryochemical potential 
Instead of using the hadronization temperature as scale, the vacuum pressure can be deter- 



mined in terms of the gluon condensate (G 2 ) [Ei|, 26], 



B 



33 - 2N f 
3847T 2 



(Gl v ). 



(5.5) 



Numerical studies [25] lead for two flavours to B 1 / 4 ~ 0.21 —0.25 GeV; from this we get both 
T Q ~ 0. 150 - 0. 175 GeV and Hq ~ 1.3 - 1.6 GeV. 

Clearly this approach is overly simplistic, since it reduces the effect of the additional quarks to 
only their Fermi repulsion. A more general way of addressing the problem would be to introduce 
an effective /I -dependence of the string tension. The presence of further quarks will lead to a 
screening-like reduction of the force between a given QQ pair, and hence a screened string tension, 
with a jit -dependent screening mass as obtained in finite density lattice studies, might provide a 
more realistic approach to the jj. dependence of the Hawking-Unruh hadronization temperature. 
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6. Angular Momentum and Non-Central Collisions 

The dependence of Hawking radiation on the angular momentum of the emitting system intro- 
duces a particularly interesting aspect for the "white hole evaporation" we have been considering. 
Consider a nucleus-nucleus collision at non-zero impact parameter b. If the interaction is of col- 
lective nature, the resulting system can have some angular momentum orthogonal to the reaction 



plane (see Fig. [llj). For central collisions, this will not be the case, nor for extremely peripheral 
ones, where one expects essentially just individual nucleon-nucleon collisions without collective 
effects. 




Figure 11: Rotating interaction region in non-central AA collision 



If it possible to consider a kinematic region in which the interacting system does have an 
overall spin, then the resulting Hawking radiation temperature should be correspondingly reduced, 



as seen in eq. ( |4.8| ). The effect is not so easily quantified, but simply a reduction of the hadronization 
temperature for non-central collisions would quite indicative. Such a reduction could appear only in 
the temperature determined by the relative abundances, since, as we shall see shortly, the transverse 
momentum spectra should show modifications due to the role of the ergosphere. The hadronization 
temperature determined from a resonance gas study could thus start for central collisions with the 



value (3.19) or, at finite fx with the corresponding Th(h), then decrease with the onset of white 
hole rotation for non-central collisions, and finally increase again as collective nuclear effects go 
away and we recover elementary nucleon-nucleon collisions. 

We next turn to the momentum spectrum of the Hawking radiation emitted from a rotating 
white hole. As discussed in section 4, such radiation will exhibit an azimuthal asymmetry due 
to the presence of the ergosphere, which by its rotation will affect the momentum spectrum of 
any passing object. At the event horizon, the momentum of all radiation is determined by the 
corresponding Hawking temperature dO|); but the passage of the ergosphere adds rotational motion 
to the emerging radiation and hence increases its momentum. As a result, only radiation emitted 
directly along the polar axis will have momenta as specified by the Hawking temperature; with 



increasing latitude 6 (see Fig. |12|)a, the rotation will increase the radiation momentum up to a 
maximum value in the equatorial plane. 

Hawking radiation from a rotating source thus leads for nuclear collisions quite naturally to 
what in hydrodynamic studies is denoted as elliptic flow. It is interesting to note that both scenarios 
involve collective effects: while in hydrodynamics, it is assumed that non-central collisions lead to 
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Figure 12: Transverse plane view of a non-central AA collision 

an azimuthally anisotropic pressure gradient, we have here assumed that such collisions lead to an 
overall angular momentum of the emitting system. 

7. Stochastic vs. Kinetic Thermalization 

In statistical mechanics, a basic topic is the evolution of a system of many degrees of freedom 
from non-equilibrium to equilibrium. Starting from a non-equilibrium initial state of low entropy, 
the system is assumed to evolve as a function of time through collisions to a time-independent 
equlibrium state of maximum entropy. In other words, the system loses the information about its 
initial state through a sequence of collisions and thus becomes thermalized. In this sense, thermal- 
ization in heavy ion collisions was studied as the transition from an initial state of two colliding 
beams of "parallel" partons to a final state in which these partons have locally isostropic distribu- 
tions. This "kinetic" thermalization requires a sufficient density of constituents, sufficiently large 
interaction cross sections, and a certain amount of time. 

From such a point of view, the observation of thermal hadron production in high energy col- 
lisions in general, in e + e~ and pp interactions as well as in heavy ion collisions, is a puzzle: how 
could these systems ever "have reached" thermalization? Already Hagedorn had therefore con- 



cluded that the emitted hadrons were "born in equilibrium" [27]. 

Hawking radiation provides a stochastic rather than kinetic approach to equilibrium. The bar- 
rier to information transfer provided by the event horizon requires that the resulting radiation states 
excited from the vacuum are distributed according to maximum entropy, with a temperature deter- 
mined by the strength of the "confining" field. The ensemble of all produced hadrons, averaged 
over all events, then leads to the same equilibrium distribution as obtained in hadronic matter by 
kinetic equilibration. In the case of a very high energy collision with a high average multiplicity 
already one event can provide such equilibrium; because of the interruption of information transfer 
at each of the successive quantum colour horizons, there is no phase relation between two succes- 
sive production steps in a given event. The destruction of memory, which in kinetic equilibration 
is achieved through sufficiently many successive collisions, is here automatically provided by the 
tunnelling process. 

So the thermal hadronic final state in high energy collisions is not reached through a kinetic 
process; it is rather provided by successively throwing dice. 
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8. Conclusions 

We have shown that quantum tunnelling through the colour confinement horizon leads to ther- 
mal hadron production in the form of Hawking radiation. In particular, this implies: 

• The radiation temperature 7# is determined by the transverse extension of the colour flux 
tube, giving 



in terms of the string tension a. 

• The multiplicity v(s) of the produced hadrons is determined by the increase of the flux tube 
thickness with string length, leading to 



where yfs denotes the collision energy. 

• The temperature of Hawking radiation can in general depend on the charge and the angular 
momentum of the emitting system. The former provides the dependence of the hadronization 
temperature on baryon density. Given sufficiently much collective collision behaviour, the 
latter could cause elliptic flow and a centrality dependence of 7# in AA collisions. 

• In statistical QCD, thermal equilibrium is reached kinetically from an initial non-equilibrium 
state, with memory destruction through successive interactions of the constituents. In high 
energy collisions, tunnelling prohibits information transfer and hence leads to stochastic pro- 
duction, so that we have a thermal distribution from the outset. 
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